Stable fractional flux vortices in mesoscopic superconductors 
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Conventional superconductors have vortices carrying integer multiples of magnetic flux quantum 
while unconventional ones, with p- or d-wave order parameter, allow half-integer fluxes. Here we 
show that mesoscopic size effects stabilize fractional flux vortices in the thermodynamical ground 
state of s-wave two-gap superconductors. The value of these fluxes can be an arbitrary fraction of 
flux quantum and can be measured directly from distributions of magnetic fields on the samples. 



PACS numbers: 03.75.Mn, 74.25.-q, 74.78.Na 



The quantization of magnetic fluxes associated to vor- 
tices is one of the most direct demonstrations of the ex- 
istence of macroscopic quantum coherence in supercon- 
ductors. The phase of superconducting order parameter 
(OP) can change only by multiples of 2ir along any closed 
contour in ordinary superconductors, which implies that 
the magnetic flux threading the vortex is quantized in 
units of $o = hc/2\e\. In contrast, the spin degree of 
freedom in triplet OP topologically allows half-quantum 
vortices^. However the energy of these vortices diverges 
in bulk samples due to unscreened spin currents^. The- 
oretical arguments were given for their stabilization by 
formation of pairs or by mesoscopic size effects, where 
the divergent part of the energy is removed- However 
the question whether such phases could be realized as 
true thermodynamical ground state remains open. Half- 
quantum vortices have been also predicted and actually 
observed in heterostructures by combining d-wave super- 
conductors with superconductors of other symmetries, 
which resulted in 7r-Josephson junctions^&i. 

Another type of fractional flux vortex was sug- 
gested in the context of multicomponent/multiband 
superconductors^'-. Multicomponcnt models gained re- 
cently much attention in connection with the discov- 
ery of two-gap superconductivity in MgBy^ and iron 
pnictidesii, where they are expected to exhibit a plethora 
of new propertie o 12 i 13 ' 14 i 15 i 16 i 17 , and also appear in the 
description of the superfluidity of liquid metallic hydro- 
gen or deuterium^ and Bose-Einstein condensates of di- 
lute atomic gasesi£. Fractional vortices in these systems 
are quite different from those arising in unconventional 
superconductors since they are neither bound to inter- 
faces, like the n-j unction vortices, nor are they related 
to a special symmetry of OP. However, as in the case 
of unconventional superconductors, such vortex phases 
are energetically forbidden in the bulk^. This restric- 
tion is in full accord with the fact that no fractional flux 
vortices have been observed to date in all investigated 
multicomponent superconductors. Contrary to that, we 
prove here that fractional flux vortices exist in mesoscopic 
two-component superconductors as thcrmodynamically 
stable phases in a broad range of temperatures, applied 
magnetic fields and superconducting parameters. More- 
over, our calculations show that fluxes associated to these 



vortices can take arbitrary fractions of $0 and should be 
observable with existing measuring techniques of mag- 
netic field distributions. 

For arbitrary (non-periodic) vortex patterns, the flux 
associated to each vortex threads the domain confined by 
a contour T which obeys the relation 
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where ^(r) is the Ginsburg-Landau (GL) OP and j(r) is 
the current density. Such definition of T automatically 
implies the quantization of flux associated to a vortex 
in conventional superconductors 1 —. In mesoscopic super- 
conductors, the quantization of the flux associated to a 
vortex is only defined when the contour T lies entirely 
within the sample. When this condition is not fulfilled 20 , 
the fractionalization of this flux arises in conventional 
superconductor s 21 ! 22 as a simple consequence of the fact 
that it has not entered yet entirely in the sample. Note 
that even if the fluxes of individual vortices are quantized, 
the total flux threading the mesoscopic sample will be 
not. This is because a non-quantized flux passes through 
the domain between the contour(s) T and the boundary 
of the sample, where the current flow is diamagnetic. To 
be distinguished from these situations, here we consider 
the case of fractionalization of the flux associated to a 
vortex when it has entered completely in the supercon- 
ductor, i.e. the contour T lies entirely within the sample. 

The current density for the 
superconductor is given by the 
j = -2efi[Im(**V*i)/mi 
4(e 2 /c)(|* 1 | 2 /m 1 + |* 2 | 2 /m 2 )A 
the OPs in the two condensates, 



two-component 
GL expression 

Im(*2V*2)A"2] - 

where ^/i and 'J 2 are 
2e is the charge of the 



Cooper pair and A is the vector potential corresponding 
to the magnetic induction b(r). Choosing a contour T 
satisfying the relation 
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we obtain for the magnetic flux threading the region con- 
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fined by T: 
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For ^(r) oc ^i(r) Eq. ([3|) reduces to the familiar flux 
quantization, while in the case of a composite vortex 
with winding numbers L\ and L 2 in the two conden- 
sates, respectively, we obtain in the London limit*: <&r = 
(Lil^lVmi + L 2 \^ 2 \ 2 /m 2 )/(\^> 1 \ 2 /m 1 + l*2| 2 /w 2 )$o. 
The last expression shows that the flux can become frac- 
tional if L\ ^ L 2l which, however, is never the case in 
bulk superconductors as mentioned above. 

Consider a long superconducting cylinder of radius 
R ~ £(0) in external homogeneous field H applied 
along its axis, described by a two-band GL free energy 
functional 2 ^: 
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-J [^( b ~ H ) 2 -7(*i*2+M5)]dV; (4) 



where the single-band functional F n = J (a n \^ n \ 2 + 

/3„|* n | 4 /2 + \n^ n \ 2 /2m n )dv with n = -mv + 

^ A are supplemented with the magnetic energy 
and the Josephson-type coupling of the two conden- 
sates. For long cylinders the problem reduces to 
two dimensions and the solutions are presented in 
t he form : ¥ a = l/ft El u L (r/R)e iLe , * 2 = 

^\ai\//3iJ2L v L(r/R)e tLd . The stable state is found by 
full numerical minimization 24 of the GL functional with 
respect to ul, Vl and A in each point of space under 
the boundary conditions n\&i| n = 0, IM^U = 0, and 
b| t = H at r = R. 

We will assume that the first band is dominant, i.e. 
a 2 becomes negative at a temperature lower than a,\. 
The coherence length is related to temperature via £ 2 = 
— h 2 /2mia 1 and a 2 /|ai| = x 20 (^i/R) 2 — ft/ft where x 20 
is a constant that quantifies the difference of intrinsic crit- 
ical temperatures between the two bands. It is then con- 
venient to parametrize the phase diagram with (i?/£i) 2 
and $ = irR 2 \H\. Fig. 1 shows the calculated phase dia- 
gram for m\ — m 2 and f3\ = /3 2 at different values of GL 
parameter n (= n\ = k 2 = $omij3^ 2 /(2ir) 3 / 2 h 2 ) and 
interband coupling x-y — j/(h 2 /2miR 2 ). At high tem- 
perature (close to the T C (H) boundary line) the OPs are 
characterized by a giant vortex in the center of the cylin- 
der (Li = L2) much as in the case of conventional super- 
conducting cylinder s 25 ' 26 . We can see that a cylinder of 
mesoscopic radius is an effectively type-II superconduc- 
tor (vortices penetrate it) even for k < 1/V2 (Fig. 1, 
bottom), like in superconducting thin films^. However 
for values n < 0.35 -j- 0.05 the whole superconducting 
domain on the phase diagram is occupied by the Meiss- 
ner phase (Vo). The shaded closed regions (Vl .l+i) m 
Fig. 1 correspond to domains with fractional flux vor- 
tices as stable thermodynamical phases. Fig. 2 shows 
indeed that the phases with fractional flux vortices are 
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FIG. 1: (color online) Phase diagram of the two-component 
superconducting cylinder for k = 10 (top), 3 (middle) and 0.5 
(bottom) (x2a = 12, 20 and 0.5, respectively) with x-y =0.01. 
It is divided into domains of superconducting states with no 
vortex (Vo), with a central giant vortex of winding number 
L in each condensate (CVi), and with L n separated vortices 
in the n-th condensate (Vl 1 ,l 2 )- Dashed and dotted lines 
delineate the domains where fractional flux vortices exist as 
stable phases for x-, =0.05 and 0.1, respectively. The stars 
on the plots denote points at which current distributions are 
shown in Fig. 3. 



the lowest in energy in the corresponding domains. As 
seen in Fig. 1, these domains shrink with increasing 7 
and decreasing k. 

In order to calculate the flux associated with vortices 
$r, we find first the contour T using the calculated cur- 
rent density distributions (Fig. 3 c-f) and Eq. ([2]). To 
this end we draw a closed line through the points where 
j(r) has zero tangential component. Then <&r is found 
either by applying Eq. ([3]) or by integrating the mag- 
netic induction threading the contour. The resulting $r 
is shown in the lower panel of Fig. 2 as function of ap- 
plied flux. We can see from the figure that in the domains 
where the $r is fractional, it almost does not varies with 
the field, looking like a step of intermediate (fractional) 
height. The height of these steps evolves monotonically 
from zero at some higher temperature till some fractional 
value at T=0. The overall picture is shown in Fig. 4. It 
shows, in particular, the existence of finite domains on 
the phase diagram where the fractional flux varies mono- 
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FIG. 2: (color online) Top: the free energy of the two- 
component superconducting cylinder for different symmetries 
of the OP for k = 10 and R 2 /£,l = 25. Vertical dot lines 
delineate the domains with fractional flux vortices in the sta- 
ble thermodynamic phase. Bottom: The flux associated to 
vortices ('I'r) in the stable thermodynamic phase. 



tonically in a broad range of values. Thus continuous 
changes of fractional flux could be detected via the mea- 
surements of the field distributions at different tempera- 
tures. 

To understand the reason for the fractionalization of 
$r we investigate the OPs and \&2 in different regions 
of the phase diagram (Figs. 3a, b). At a high tempera- 
ture, (dashed line) the two OPs correspond to the giant 
vortex with the same vorticity (L — 1) and moreover are 
proportional to each other. Then according to Eq. ([3]) we 
obtain $r/^o = 1, i-e. the vortex flux is quantized. For 
temperature low enough the giant vortex in the second 
band mostly corresponds to L = 2, while the active band 
still contains the same giant vortex with L = 1. Given 
this difference in the two OPs, Eq. ([3]) will result in a 
fractional value of <i>r. The current distribution for this 
state is shown in Fig. 3d. 

The Josephson-type interaction between two conden- 
sates favors maximal overlap of ^i(r) and ^(r)- The 
two condensates also interact indirectly through the 
magnetic induction. When expressed with normal- 
ized quantities, it can be shown to the first order in 

K~ 2 that F = E n ^n[* n ,A ] - 7/(*l*2 + *1*2) - 

K ~ 2 II J'( r )j( r ')/I r — r 'li where F n [^ n , A ] is the func- 
tional of the n-th condensate in the unperturbed applied 
field with the vector potential Ao. The last term in the 
normalized functional favors the maximal overlap of the 
current distributions in the two condensates, i.e. the 
maximal similarity of the two OPs. This explains why 
the domains of the existence of fractional flux vortices 
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FIG. 3: (color online) OP components ul (a) and vl (b) for 
R 2 /£,! = 24 (dashed line) and 27 (solid line), corresponding to 
quantized and fractional vortex phases, (c) Radial variation 
of the magnetic induction in the corresponding phases, (d) 
Current density distributions and streamlines for values of 
-R 2 /£i indicated by stars on Fig. 1. The position of vortices 
in the two condensates are denoted by crosses for the first 
band and by squares for the second band. The contour V is 
shown by the white dashed line. 



shrink with the increase of 7 (x 7 ) and the decrease of k 
(Fig. 1). 

For MgB2, the most investigated two-band supercon- 
ductor, the Josephson coupling is moderately weak: the 
smallest estimate is 7 > 0.3|ai which yields x 7 > 
0.3(i?/£i) 2 . Our numerical minimizations have not found 
any flux fractionalization for x 1 > 0.2 or (i?/£i) 2 < 3. 
These limitations mean that 7 in MgB2 is at least five 
times too large. However there is no reason why 7 
cannot be lower in other multiband superconductors. 
Among the promising candidates we can cite the heavy 
fcrmion compound UNi2Ai3— , the quaternary borocar- 
bide YN^P^C^ii or the ternary-iron silicide Lu2Fc3Si5 31 
where stronger deviations from the predictions for a sin- 
gle gap suggest that their interband couplings may be 
weaker than in MgB2. 

The flux fractionalization can be understood as re- 
sulting from the combination of two weakly interacting 
condensates which, when uncoupled, individually prefer 
states of different vorticities for the same applied mag- 
netic field. So this mechanism is not restricted to the 
particular sample symmetry studied here and should be 
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present in other geometries too. For example the unlock- 
ing of the condensate vorticities is not specific to long 
cylinders since it exists as well in the opposite limit of 
thin disks. We have furthermore checked that it can also 
occur in cylinders with different central holes^. 

We notice from Fig. 3 that the contour T is circular 
for central vortex phases and is still very close to a circu- 
lar shape for OPs characterized by low vorticity, despite 
the fact that the vortices are deconfined (vortices in the 
two condensates are spatially separated). The distribu- 
tion of the magnetic induction for these states (Fig. 3c) 
shows a characteristic depletion at a radius which coin- 
cides exactly with the position of the contour T. There- 
fore for experimental determination of <£>r in these states 
it would be sufficient to integrate the measured distribu- 
tion of b(r) in the region confined by this depletion. The 
amplitude of the variation of magnetic induction is ca 5 
% for k = 1 and ca 20 % for k — 3, which is larger than 
the accuracy of the current measuring techniques. For 
phases of higher vorticity or samples of non-cylindrical 
geometry the contour T is not circular anymore. In such 
cases the distribution of currents should be obtained first 
from the measured distribution of magnetic induction on 
the top surface using the relation j(r) = (c/4ir)V x b(r). 



Then the contour T can be found as the line on which 
the current has no tangential component. 

In conclusion, we have proven by exact calculations 
within GL theory that fractional vortices, carrying ar- 
bitrary fraction of flux quantum (Fig. 4), exist in meso- 
scopic two-component superconductors as stable thermo- 
dynamic phases in a broad range of superconducting pa- 
rameters, temperature and applied magnetic field. The 
origin of stabilization of these vortex phases is the differ- 
ent effect of the confinement on the two condensates ex- 
erted by the boundary of the sample, i.e. is a pure meso- 
scopic effect. This effect is counteracted by the Joscphson 
coupling between the two condensates and the magnetic 
field induced by the screening currents, both favoring 
similar distributions of the OPs. The obtained fractional 
flux vortices should be observable with available experi- 
mental techniques for the investigation of magnetic field 
distribution. 
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